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Abstract: We present forecast errors on a wide range of cosmological parameters obtained 
from a photometric cluster catalogue of a future wide-field Euclid-like survey. We focus 
in particular on the total neutrino mass as constrained by a combination of the galaxy 
cluster number counts and correlation function. For the latter we consider only the shape 
information and the Baryon Acoustic Oscillations (BAO), while marginalising over the 
spectral amplitude and the redshift space distortions. In addition to the cosmological 
parameters of the standard ACDM+z^ model we also consider a non-vanishing curvature, 
and two parameters describing a redshift evolution for the dark energy equation of state. 
For completeness, we also marginalise over a set of "nuisance" parameters, representing the 
uncertainties on the cluster mass determination. Wc find that combining cluster counts 
with power spectrum information greatly improves the constraining power of each probe 
taken individually, with errors on cosmological parameters being reduced by up to an order 
of magnitude. In particular, the best improvements are for the parameters defining the 
dynamical evolution of dark energy, where cluster counts break degeneracies. Moreover, 
the resulting error on neutrino mass is at the level of cr(Mj,) ~ 0.9 eV, comparable with that 
derived from present Lya forest measurements and Cosmic Microwave background (CMB) 
data in the framework of a non-flat Universe. Further adopting Planck priors and reducing 
the number of free parameters to a ACDM-|-z/ cosmology allows to place constraints on the 
total neutrino mass of a{My) ~ 0.08 eV, close to the lower bound enforced by neutrino 
oscillation experiments. Finally, in the optimistic case where uncertainties in the calibration 
of the mass-observable relation were so small to be neglected, the combination of Planck 
priors with cluster counts and power spectrum would constrain the total neutrino mass 
down to (t{M,j) ~ 0.034 eV, i.e. the minimum neutrino mass predicted by oscillation 
experiments would be detected in a ACDM framework. We thus show that galaxy clusters 
from future wide galaxy surveys will be an excellent tool for studying cosmology and 
fundamental physics. 
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1. Introduction 

It is now established from solar, atmospheric, reactor and accelerator experiments that 
neutrinos have non-zero mass, and that a lower limit on the total neutrino mass is given 
by Ml, = Y^rrii, ~ 0.05 eV [1, 2], where rrii, is the mass of a single neutrino species. On 
the other hand the individual masses are still unknown. Since neutrino mass affects the 
evolution of the Universe in several observable ways, its measurements can be obtained 
from different cosmological probes, such as observations of the CMB, galaxy clustering, 
Lya forest, and weak lensing data [3, 4, 5, 6]. 

In particular, a thermal neutrino relic component in the Universe impacts both the ex- 
pansion history and the growth of cosmic structures. Neutrinos with mass < 0.6 eV become 
non-relativistic after the epoch of recombination probed by the CMB, and this mechanism 
allows massive neutrinos to alter the matter-radiation equality for a fixed Q,„ih^- Massive 
neutrinos act as non-relativistic particles on scales k > knr = 0.018(m,//leV)^/^rim^/i/Mpc, 
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where /cnr is the w8,ve-nuiiiber corresponding to the Hubble horizon size a,t the epoch z^i^ 
when the given neutrino species becomes non-relativistic, Q^n is the matter energy density 
and h = i7o/100kms~^Mpc~^. The large velocity dispersion of non-relativistic neutri- 
nos suppresses the formation of neutrino perturbations in a way that depends on rrii, and 
redshift z, leaving an imprint on the matter power spectrum for scales k > kfs{z) = 
0.82H{z)/Hq/{1 + z)^(mjy/leV)/i/Mpc [1, 7], where neutrinos cannot cluster and do not 
contribute to the gravitational potential wells produced by cold dark matter and baryons. 
This modifies the shape of the matter power spectrum and the correlation function on 
these scales [8, 9, 10, 11, 12, 13, 14]. 

Massive neutrinos affect also the CMB statistics. WMAP7 alone constrains < 1.3 
eV [15] and data from the ACT^ and SPT^ experiments constrain M^, < 0.948 cV in the 
framework of a non-flat cosmology [16]. FTirthcrmore, thanks to the improved sensitivity 
to polarisation and to the angular power spectrum damping tail, forecasts for the Planck 
satellite alone give a 1-a error on the total neutrino mass of ~ 0.2 — 0.4 eV, depending 
on the assumed cosmological model and fiducial neutrino mass [17, 18]. Moreover, the 
combination of present data-sets from CMB and large-scale structure (LSS) yields an upper 
limit of Mjy < 0.3 eV [19, 20, 21, 22, 23]. A further robust constraint on neutrino masses 
has been obtained using the Sloan Digital Sky Survey flux power spectrum alone, finding 
an upper limit of < 0.9 cV {2a C. L.) [24]. However, the tightest constraints to date 
in terms of a 2a upper limit on the neutrino masses have been obtained by combining the 
Sloan Digital Sky Survey flux power from the Lya forest with CMB and galaxy clustering 
data. These constraints are partly driven by a discrepancy in the measured erg between 
Lya forest and CMB, and result in Sm,y < 0.17 eV [25]. Somewhat less constraining 
bounds have been obtained by [26], while for forecasts on future CMB and Lya forest joint 
constraints we refer to [27]. 

Moreover, the forecast sensitivity of future LSS experiments, when combined with 
Planck CMB priors, indicates that observations should soon be able to detect signatures 
of the cosmic neutrino background and measure the neutrino mass even in the case of 
the minimum mass = 0.05 eV [18, 28, 29, 30, 31, 32]. Furthermore, these future 
surveys have been planned to measure with high accuracy the so-called "dark energy" 
equation of state. In fact, there is now strong evidence that the current energy density 
of the Universe is dominated by dark energy with an equation of state it; ~ — 1, which 
is causing accelerated expansion. Accordingly, in this work we consider, in addition to 
massive neutrinos, the effect on LSS of a homogeneous dark energy component with a 
general time- varying equation of state w{z). 

In this respect, Ref. [32] has shown that future spectroscopic galaxy surveys, such as 
the ESA selected mission Euclid^ [33], will be capable of estimating the neutrino mass 
scale independently of flatness assumptions and dark energy parametrisation, if the total 
neutrino mass M,y is > 0.1 eV. On the other hand, if M,y is < 0.1 eV, the sum of neutrino 
masses, and in particular the minimum neutrino mass required by neutrino oscillations, 

^http://wwwphy.princeton.edu/act/ 

^http:/ /pole. uchicago.edu/ 
^http://www. euclid-ec.org/ 
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can be measured in the context of a ACDM model. For further discussion on neutrino 
mass constraints from different probes see e.g. [2, 3] and references therein. 

In this picture, httle attention has been dedicated to the constraining potential of 
galaxy clusters. Historically these systems have played a fundamental role in cosmology. 
For instance, the correlation function of galaxy clusters provided the very first hint toward 
the low-density Universe that is today commonly accepted as the standard cosmological 
model [34, 35]. Galaxy cluster number counts at high-redshift later provided yet another 
early evidence for the matter density parameter being substantially smaller than unity 
[36, 37]. Galaxy clusters trace the large-scale matter distribution much as galaxies, with 
the difference that the former have a substantially larger bias than the latter, and hence 
their correlation is a factor of a few higher. At the same time, clusters are relatively rare 
objects, therefore their abundance is highly sensitive to features in the primordial matter 
power spectrum. These occurrences make clusters in principle valuable tools for constrain- 
ing the details of the cosmological model, and in particular the neutrino masses. However 
their low spatial number density is also a drawback, in that it increases substantially the 
shot noise and the Poisson noise as compared to galaxies. Their effective cosmological 
power hence depends significantly on the chosen selection function and on the precision 
with which cluster masses can be estimated. In this work we study how the combination 
of cluster number counts with the shape of the cluster power spectrum, including infor- 
mation from the BAO, can constrain the total neutrino mass. We focus on the cluster 
catalogue that will be produced by the photometric part of Euclid, and provide constraints 
on the parameter set of an extended cosmological model which includes curvature, massive 
neutrinos, and a time- varying dark energy component. Since deviations from the standard 
model in the form of extra neutrino species are still uncertain [38, 39], in this paper we focus 
on standard neutrino families only and analyse the constraining power that future galaxy 
cluster catalogues, photometrically selected from galaxy surveys of Euclid-type, have on 
the total neutrino mass. 

The rest of the paper is organised as follows. In § 2 we present the galaxy cluster 
modelling adopted in this work and the modifications due to uncertainties in the calibration 
of the mass-observable relation. In § 3 we describe the adopted fiducial model and the 
effect of massive neutrinos on the matter power spectrum. In § 4 we review our forecasting 
approach as applied to cluster counts and power spectrum, and specify characteristics of 
the galaxy cluster survey analysed in this work. In § 5 we present our results on the forecast 
neutrino mass errors, and finally in § 6 we draw our conclusions. 

2. Clusters: effective bieis and mass function 

In the cosmological analysis performed in this paper we combine number counts and the 
power spectrum of galaxy clusters. These are identified as overdensities of galaxies pho- 
tometrically selected with a Euclid-type survey. In order to properly define the relevant 
cluster catalogue we adopt the minimum mass Minin(z) provided by the Euclid Red-Book 
[33] for objects identified as having a S/N ratio larger than 5. Let us assume a perfect 
knowledge of the mass of each cluster in the catalogue for the time being. This assump- 
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tion will be relaxed in § 2.1 below. The resulting (full-sky equivalent) cluster redshift 
distribution can then be computed as 



dN{z) _ dV{z) 



9{z) , 



(2.1) 



dz dz 



where dV{z)/dz is the cosmic volume per unit redshift, while 




Mniin{2) 



dM n{M, z) . 



(2.2) 



For the mass function n(M, z) of dark matter halos we adopt the prescription given by 
[40], which is based on an approximated ellipsoidal collapse model and has been shown to 
agree well with the results of numerical cosmological simulations. Although derived in a 
ACDM context, this prescription has been shown to model n(M, z) accurately even in the 
presence of dynamic dark energy [41, 42] and massive neutrinos [43, 14, 44], and holds, in 
particular, for the fiducial values of neutrino mass and dark energy parameters considered 
in this work (see §3). 

As for the power spectrum of galaxy clusters, we simply assume it to be a biased 
version of the linear dark matter power spectrum. Nonlinear effects become important on 
scales k > 0.1 h Mpc^^, which at low redshift have been disregarded in our Fisher matrix 
analysis, as wc shall explain in Appendix A. By assuming again a perfect knowledge 
of the cluster masses in the Euclid catalogue, we can compute the power spectrum as 
P{k,z) = b1{z) Pi^{k,z), where Pi{k,z) is the linear dark matter power spectrum, while 
be{z) is the effective bias of clusters in the catalogue, defined as 



In Eq. (2.3), the function b{M,z) represents the bias of dark matter halos, for which we 
adopt the semi-analytic prescription of [45] (see [43] for an analysis of this prescription 
against N-body simulations which include a massive neutrino component). We compute 
the linear dark matter power spectrum Phik, z) with the publicly available software package 
CAMB [46], which takes correctly into account the effect of massive neutrinos. 

2.1 Nuisance parameters 

The discussion presented in § 2 assumes a perfect knowledge of the true mass of clusters 
in our catalogue. This is a very strong assumption, in that reliable mass estimates can be 
obtained only with an extensive multi- wavelength follow-up and only for the most massive 
objects. A less expensive alternative, particularly suitable for large cluster catalogues, is 
to adopt a scaling relation between the observable quantity at hand (cluster richness in 
this case) and the true underlying mass. Scaling relations between cluster properties are 
however not perfect one-to-one associations; they include a scatter and, in some circum- 
stances, systematic biases. In order to take properly into account the uncertainties that 
a scaling relation introduces in the knowledge of the true underlying mass, we treat the 




+ 0O 



(2.3) 
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scatter and systematic biases as "nuisance" parameters, and marginalise over them in the 
Fisher matrix analysis for cluster counts in § 4.2. 

Given this, let p{Mo\M) be the probability that for a given cluster of intrinsic mass M 
we infer the mass Mq through the scaling relation [47] . The redshift distribution of clusters 
given by Eq. (2.2) above then gets modified according to 



r+oo r+oo 

g{z)= dMo j dM n{M,z) p{Mo\M). (2.4) 

Jm^^^{z) Jo 



By assuming a lognormal scatter around the nominal scaling relation with dispersion ainM, 
the probability can be written as 

P(^o\M) = -— i exp i-x^M,)] , (2.5) 

MoV^TTCJlnM 

where ^ 

x(Mo) = [In(Mo) -Bm- ln(M)] . (2.6) 

V 2o-lnM 

Evidently, the parameter Bm represents the fractional value of a systematic bias in the 
scaling relation. It easily follows that the function g{z) can then be simplified to 

1 r+oo 

g{z) = - dM niM, z) erfc [x{M^in{z))] , (2.7) 
^ Jo 

where erfc(x) is the complementary error function [48]. 

Likewise, the effective bias presented in Eq. (2.3) now takes the more general form 

1 f+oo 

b,{z) = — - / dM n{M, z) b{M, z) erfc [x(M„i„(z))] . (2.8) 

^g{z) Jo 

It is straightforward to verify that, in the ideal case in which the systematic bias Bm 
and the scatter (XinM ~^ 0, the expressions for the cluster redshift distribution and the 
effective bias reduce to their simplest form, presented at the beginning of §2. In what 
follows we adopt these updated forms for the cluster redshift distribution and the effective 
bias. 

Moreover, following [47], we assume the following redshift parametrisation for the halo 
mass bias and variance: 

\nBM{z) = A + Bln{l + z) 

^lnM(^) = ^lnM,0 " 1 + (1 + zf^ . (2.9) 

In this way, we have four nuisance parameters, namely A, B, (TinM,o aiid /3. For the reference 
model, we assume In Bm{z) to be zero at 2; = (i.e. A = 0), with no evolution (i.e. B = 0), 
but leaving A and B as free parameters with Gaussian priors of o-{A) = a{B) = 0.25 [49], 
which is not overly restrictive for Euclid [33]. In Eq. (2.9), we have also allowed for a 
power-law redshift evolution for the variance crf^j^{z) of the lognormal intrinsic scatter 
around the nominal scaling relation. Ref. [50] estimates that cri^M,o = 0.2 and we choose 
P = 0.125. This means that the scatter will grow to a value of 0.6 at a redshift of 2: = 2. In 
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the count Fisher matrix analysis we self-calibrate for these scatter variables with Gaussian 
priors of cr{ainMfl) = 0.1 [51] and cr{(3) = 0.1, which should be conservative estimates, 
although these quantities are unconstrained by present data [33]. It is expected that Stage 
III surveys, e.g. DES^, will shed light on this quantity beyond redshift z = 1. 

As mentioned above, for the minimum cluster mass as a function of redshift Mminiz), 
we adopt the prescription detailed in the Euclid Red-Book [33] for the photometric cluster 
catalogue. The redshift evolution of the cluster number counts is computed in bins of width 
Az = 0.1 between z = 0.2 and z = 2, and we integrate the cluster distribution over mass, 
considering number counts in redshift-space only. Concerning redshift errors that enter the 
observed galaxy cluster power spectrum, better described in Appendix A, we assume that 
photometrically selected clusters will be spectroscopically confirmed up to 2; = 1. 

3. Cosmological model with massive neutrinos 

According to the latest observations (e.g. [15, 52] and refs. therein), we assume the 
following fiducial cosmological model at the present epoch: = 0.271, = 1 — 
h = 0.703, As = 2.525X 10-9, = 0.045, = 0.966, wq = -0.95, Wa = 0,M^ = = 
0.05eV^. This corresponds to cg = 0.8 for the fiducial cosmology. We consider neither 
primordial gravitational waves nor a scale dependent component of the scalar spectral 
index, and assume the matter energy density to include the baryon and neutrino 
contributions when neutrinos are non-relativistic, so that = Clc + i^b + ^u, where 
= M^/{93.Uh'^eV) [53]. 

In our fiducial cosmology, the dark energy is described by a cosmic fluid with an 
equation of state Wddz) = Pde{z) / pde{z), where and pde are the pressure and energy 
density of the dark energy fluid, respectively. The redshift dependent dark energy density 
is then 



which we normalise so that = at the present epoch. Finally, to compute our 
forecasts on dark energy parameters, we adopt the widely used linear dark energy equation 
of state Wdeia) = wq + (1 — a)wa [54, 55], where a = 1/(1 + 2;) is the scale factor normalised 
to unity at present. When the dark energy equation of state is a function of redshift, as 
we assume in the present work, the constraints on the the sum of neutrino masses can 
degrade significantly, since dark energy and massive neutrinos both affect the growth rate 
of structures [56]. However, as we will show in § 5, the combination of CMB and galaxy 
cluster data reduces or even breaks these degeneracies. 

As we have mentioned in § 1, massive neutrinos suppress the matter power spectrum 
on wave- numbers above the free-streaming scale /cfg, while on very large scales they behave 
as ordinary cold dark matter. The suppression effect is encapsulated in the parameter f^^ = 

''http:/ / www.darkenergysurvey.org/ 

®Here we assume an effective number of neutrino species iVefj = 3.04 and, since neutrino oscillation 
experiments have shown that at least one neutrino species is heavier than 0.05 eV, we consider an inverted 
hierarchy where mi ~ m2, ms ~ 




(3.1) 
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^u/^rm which damps the source of matter density perturbations. In fact, the hnear growth 
rate is defined as fg = din (5/(iln a, where 5 represents the matter density perturbation 
5 = bpml pmi and pm and bpm the matter density and the overdensity, respectively. In the 
presence of massive neutrinos and a dark energy component the equation for the linear 
evolution of matter density perturbations on scales k S> fcfg can be written as [57, 41] 

^ = -/| - |l - ^[^^m(a) + (3^de(a) + l)J^<ie(a)]| jg + \^mia){\ - U), (3.2) 

where Qm{o.) = NQQ^a'^^ / H'^ (a) and il^deio-) = HqQ/^X (a) / H'^ (a) are the time-dependent 
density parameters of matter and dark energy, respectively. Here H{a) = Hq{VL^^^ + 
^mO~^ + i^i^-a"^ + OAX(a))^/^, where, for the parametrisation of the dark energy equation 
of state chosen above, we have X{a) = a~^(-'^+"'o+"'°)exp(3i(;a(a — 1)) [58]. 

Eq. (3.2) is a simplified description of the effect of massive neutrinos on the growth of 
structures, since in the presence of massive neutrinos the growth rate is not only redshift- 
dependent, but also sca/e-dependent. For this reason, in the present work we have com- 
puted fg using the semi-analytic formula of Ref. [57], as we explain in Appendix A. 



4. Forecast approach: combining cluster BAO and counts 

In this work we derive neutrino constraints combining measurements from a Euclid-Vike 
galaxy cluster catalogue with CMB measurements as obtained from a Planck-YikB CMB 
experiment^. To this aim, we adopt a Fisher matrix approach [59] that allows us to 
forecast cosmological parameter errors from LSS and CMB. In order to compute CMB 
priors we use the specifications of the Planck satellite and, as explained in Appendix 
B, we describe CMB temperature and polarisation power spectra using the parameter 
set 6 = {a;^,c<;6,a;i/, 100^5, log(10^°74s),n5,r}, where Og is the angular size of the sound 
horizon at last scattering, and r is the optical depth due to reionisation. After marginalising 
over the optical depth, we propagate the Planck CMB Fisher matrix F^^^ into the final 
sets of parameters q adopted in § 4.1, using the appropriate Jacobian for the involved 
parameter transformation. 

4.1 Fisher matrix for BAO-|-Pc(k)-shape 

In this section we apply, to the observed galaxy cluster power spectrum, the so-called "P(A;) 
method marginalised over growth information" , which exploits only the shape and the BAO 
positions of the cluster power spectrum, while marginalising over amplitude and redshift- 
space distortions. This method allows us to estimate, from the galaxy cluster catalogue, 
measurements of the cosmological parameters which characterise the underlying fiducial 
cosmology. Up to now, this method has been applied to the galaxy power spectrum in a 
number of works (see e.g., [32, 60, 61, 62, 63, 64, 65, 66, 67]). Here we present the main 
formulae describing this approach, referring to Appendix A for details. 

®www.rssd.esa.int/index.php?project=planck 
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Including redshift-space distortions and the geometrical effects due to the incorrect 
assumption of the reference cosmology with respect to the true one [60], the observed 
galaxy cluster power spectrum can be written as 

-Pobs(feref±, Kef\\,z) = ;^^^^f^^^'^c(^ref±, ^ref |h ^) + ^shot , (4.1) 

where Pc is given by Eq. (A. 4) and Da{z), H{z), k±, k\\ are defined in Appendix A. 

We divide the survey volume into redshift shells with size Az = 0.1, centred at redshift 
Zi, and choose the following set of parameters to describe Pobs{kref±jkreip^)' 

{H{zi),DA{zi),G{zi),l3{zi,k),Plf,„t,uJm,i^b,i^u,ns,h} , (4.2) 

where cOi, = ^^h^, cOm = ^mh^, ^6 = ^bh^- Finally, since the growth factor G{z), the 
effective bias bc{z), and the power spectrum normalisation Pq are completely degenerate, 
we introduce the quantity G{zi) = {Po)^'^bc{zi)G{zi) /G{zo) [62]. 

Under the assumptions discussed in Appendix A, the Fisher matrix associated to the 
observed galaxy cluster power spectrum can be approximated as [68, 69] 



/-^^ ainPobs(fc)91nPobs(fc) ,, .A^^ 

i-iJk^^. dPi dp, ^^^^^'^^ 2(2.)3 ' 

where the derivatives are evaluated at the parameter values Pi of the fiducial model, and 
VeS is the effective volume of the survey: 

2 



Purvey (4-4) 



JlcPcik,IJ.) + 1 

Here we have assumed that the comoving galaxy cluster number density ric is constant in 
position and given by Eq. (2.7). Due to azimuthal symmetry around the line of sight, the 
three-dimensional galaxy cluster power spectrum -Pobs(fc) depends only on k and /x, i.e. it 
is reduced to two dimensions by symmetry [60]. 

We do not include information from the amplitude G{zi) and the redshift space distor- 
tions P{zi, k), so we marginalise over these parameters and also over P^^^^- Then we project 
p = {H{zi),DA{zi),u!m,oJb,^i^,iT's,h} into the final sets q of cosmological parameters [70], 

q = {nm,^K,^b,h,Mi,,ns,wo,Wa,log{10^^As)} . (4.5) 

The transformation from one set of parameters to another is given by 

F^s = V — F.P^ ^ U 6) 

where F^^ is the survey Fisher matrix for the set of parameters q, and Pf-" is the survey 
Fisher matrix for the set of equivalent parameters p. 

The adopted full P{k) method marginalised over growth information does not give any 
constraint on Ag, since the normalisation of the cluster power spectrum is marginalised over. 
Therefore, the log(10^'^74s)-errors shown in § 5 come from combining cluster counts and 
CMB measurements. 
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4.2 Fisher matrix for cluster counts 

Following the approach of [71, 72, 73, 74, 75], the Fisher matrix for the number of clusters, 
Ni, within the i-th redshift bin and mass M > Minin(2:), can be written as 



where the sum over i runs over redshift intervals and the number of clusters expected in a 
survey having a sky coverage AO in the redshift range between Zi and Zj+i can be written 
as 

Here g{z) is given by Eq. (2.7), dV/dz is the cosmology-dependent comoving volume el- 
ement per unity redshift interval and solid angle, and finally, the cosmological parameter 
set Qa is given by 

q = {^rn,^K,^b,h,M„,ns,WQ,Wa,\og{lQ^^As),A,B,Ui^Mfi,P} ■ (4.9) 

Where not otherwise specified, we marginalise over the four nuisance parameters A, B, 
o'inM,0) a'Hd P to obtain the constraints on Qq, of Eq. (4.5) from cluster counts. It is worth 
noting that in Eq. (4.7) we have adopted the same approach of Ref. [75], i.e. we integrate 

the cluster distribution over mass and consider number counts in redshift-space only, in 
order to reduce the effects due to the model-dependence of the cluster mass inference. 

Having defined the Fisher matrices Faj3 for CMB, Pc{k), and Ni, respectively, the 1-a 
error on Qa, marginalised over the other parameters, is a{qa) = \l {F~^)aa- Furthermore, 
to quantify the level of degeneracy between the different parameters, we estimate the so- 
called correlation coefficients, given by 

r= , ^^"'^" ^ ■ (4.10) 

When the coefficient |r| = 1 the two parameters are totally degenerate, while r = means 
they are uncorrelated. 

InS 5 we shall evaluate cr^qa) and r both from galaxy cluster data, F^^, -F^, ■Fa^'^-^a^: 
and from their combination with CMB priors, F^^+F^^^, F^^+F^^^, F^^+F^^+F^^^. 

4.3 Euclid-like survey 

In this work we forecast neutrino constraints using count and power spectrum measure- 
ments of galaxy clusters photometrically selected by a future survey like the one planned 
by the ESA selected M-class Euclid galaxy survey. This survey will be able to push to- 
wards very high redshifts over a large area, thanks to its unique capabilities in the infrared 
which cannot be matched from the ground. Conservative estimates, based on simulated 
mock catalogues, indicate that Euclid will find of order ~ 6 x 10^ clusters with S/N bet- 
ter than 3, between z = 0.2 and z = 2.0, with 10*^ at z > 1 [33]. In this case, the 
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Table 1: Forecast l-cr errors for the cosmological parameters considered in the text and the 
corresponding correlations r with M^, for a Euclid-like experiment alone and in combination with 
Planck. 
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0.0849 


0.0039 


0.3854 




0.1008 


0.6507 


0.0072 


0.0083 


1.0042 


-0.5785 


0.0050 


-0.4842 


0.0143 


-0.6459 


0.0016 


0.0125 


h 


0.3955 


0.6583 


0.0562 


-0.0158 


7.3604 


0.5735 


0.0386 


0.4771 


0.0411 


-0.5120 


0.0125 


-0.0380 




2.5972 


1.0000 


0.2283 


1.0000 


28.513 


1.0000 


0.2564 


1.0000 


0.8905 


1.0000 


0.1853 


1.0000 


Us 


0.5332 


-0.2732 


0.0024 


-0.0666 


4.0365 


-0.6573 


0.0024 


-0.0486 


0.1452 


0.8252 


0.0024 


-0.0297 


Wo 


1.7820 


0.05685 


1.4648 


-0.0248 


0.4659 


0.0997 


0.2381 


-0.2154 


0.3157 


-0.2136 


0.2133 


0.0020 


Wa 


6.6788 


-0.2651 


5.4101 


-0.0044 


5.3758 


-0.4089 


0.8187 


0.0374 


1.0644 


-0.2882 


0.7218 


-0.2123 


log(iO^"-l.) 






0.0253 


0.9117 


2.2()31 


0.9057 


0.0287 


0.933G 


0.2899 


0.182G 


0.0222 


U.8878 



ACDM 





BAG 


BAG 


f CMB 


COUNTS" 


COUNTS" +CMB 


BAO+COUNTS" 


BAO+COUNTS" +CMB 




a 


r 


a 


r 


(T 


r 


(J 


r 


(7 


r 


(7 


r 




0.1131 


0.2559 


0.0038 


-0.8535 


0.0461 


0.7086 


0.0087 


-0.9673 


0.0252 


0.3455 


0.0036 


-0.8361 




0.0308 


-0.2015 


0.0010 


-0.7291 


0.6646 


0.0628 


0.0021 


-0.9401 


0.0107 


-0.4618 


0.0009 


-0.7042 


h 


0.0935 


-0.2121 


0.0071 


0.8008 


4.7404 


0.3998 


0.0162 


0.9561 


0.0310 


-0.4705 


0.0068 


0.7796 




1.7791 


1.0000 


0.0818 


1.0000 


4.2971 


1.0000 


0.1679 


1.0000 


0.5686 


1.0000 


0.0758 


1.0000 


ris 


0.3100 


0.7894 


0.0023 


0.2498 


1.4125 


-0.6499 


0.0023 


-0.0576 


0.0713 


0.5757 


0.0023 


0.2220 


log(10iM,) 






0.0140 


0.6771 


0.9446 


0.2623 


0.0225 


0.9122 


0.1021 


0.8754 


0.0138 


0.7522 



"The four nuisance parameters, A, B, <TinM,o, and /? have been marginalised over. 



cluster-based constraints on cosmological parameters will be limited by the understanding 
of the catalogue selection function, systematic errors and cluster mass determinations and 
their uncertainties. In this respect, Euclid will be able to calibrate the mass-observable 
relations and their scatter through lensing measurements. The high image quality and 
number density of sources will enable Euclid to measure masses of clusters much more 
accurately and out to higher redshifts than is possible from the ground. Moreover, the 
combination of Euclid data with other surveys, such as eROSITA^ or Planck, will enable 
the cross-calibration of non-lensing mass-observable relations, which are currently limited 
to low redshifts and small samples. For instance, Euclid will provide mass proxies via the 
stacking of clusters according to their X-ray or Sunyaev-Zel'dovich signals. [33]. Here we 
adopt the following specification for the Euclid-like galaxy cluster catalogue analysed in 
this work: area= 15 x lO^deg^, Zmin = 0.2, Zmax = 2, = 0.1, S/N> 5. 

5. Results 

In this Section we present the predicted l-cr marginalised errors and correlations for the 
cosmological parameter set, described in Eq. (4.5), considered in this work, focusing in 

http://www.mpe. mpg.de/erosita/ 
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Table 2: Comparison of the Euclid+ Planck forecast 1 a errors and the corresponding correlations 
with My, for different treatments of the nuisance parameters. 



General cosmology 





BAO+COUNTS»+CMB 


BAO+COUNTS^+CMB 




a 


r 


a 


r 




0.0106 


0.4046 


0.0057 


0.6860 




0.0039 


0.3854 


0.0031 


0.6551 




0.0016 


0.0125 


0.0008 


0.3456 


h 


0.0125 


-0.0380 


0.0061 


-0.3841 


M„ 


0.1853 


1.0000 


0.1064 


1.0000 


Us 


0.0024 


-0.0297 


0.0023 


0.0722 


Wo 


0.2133 


0.0020 


0.1477 


0.7178 


Wa 


0.7218 


-0.2123 


0.5119 


-0.7616 




0.0222 


0.8878 


0.0175 


0.8599 


ACDM 




DAo cm 


XIS" 1 i-.\lB 


BA( ) 


XI s'" 1 ( 'MB 




0" 


r 


G 


r 




0.0036 


-0.8361 


0.0008 


0.1401 




0.0009 


-0.7042 


0.0003 


0.5604 


h 


0.0068 


0.7796 


0.0021 


-0.3582 




0.0758 


1.0000 


0.0339 


1.0000 


Us 


0.0023 


0.2220 


0.0023 


0.5414 


log(10"'A.) 


0.0138 


0.7522 


0.0103 


0.9441 



"The four nuisance parameters, BMfi, cri-aMfi, a, and /3 have been marginalised over. 

''The four nuisance parameters, Bm,o, cinM,o, a, and /3 have been kept fixed to the fiducial values 

reported in the text. 

particular on the total neutrino mass M,y. As explained in § 3, we have assumed an 
inverted mass hierarchy with a fiducial M^, = 0.05 eV, which represents the minimum 
total neutrino mass predicted by neutrino oscillation experiments. In this respect, we have 
been very conservative, and consequently the forecast (T{My) errors should be viewed as 
upper bounds for neutrino mass measurements in the context of a Euclid-Wke galaxy cluster 
survey in combination with Planck. 

In Table 1 we show the marginalised errors for the cosmological parameters of Eq. (4.5) 
and the corresponding correlation coefficients with My, for the different probes and their 
combinations analysed in this work. Moreover, Figs. 1-2 represent the corresponding 2- 
parameter My-q^ joint contours at 68% C.L. Finally, in Tabic 2 we compare the cosmo- 
logical parameter errors and their correlations with My as obtained from galaxy cluster 
measurements in combination with CMB data, for the two different cases in which the 
nuisance parameters Bm,o, o'inM,0) S'l^d P are kept fixed to their fiducial values or are 
marginalised over. 

5.1 BAO+CMB 

Let us consider the different probes separately. The 2'^'^ and 3^'^ columns of the upper panel 
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Figure 1: 2-paraineter projected 68% C.L. contours in the M^-qa subspace with = 
ri/f, rib, h. The sohd black hne and the dotted blue line correspond to BAO+Pc(fc)-shape 
cluster data from a Euclid-\ike survey in combination with Planck priors for a general cosmology, 
and a ACDM Universe, respectively. The dashed red line, and the dot-dashed orange line are 
obtained with the further addition of cluster counts data, and represent the confidence contours 
obtained for the two cosmological models, respectively. 



of Table 1 show the 1-a errors and the values of the r coefficients for the Mj^-qa correlations 
obtained in a general cosmology which includes dark energy and curvature, and exploiting 
information from the combination of BAO-|-P(A:)-shape extracted from the observed galaxy 
cluster power spectrum of Eq. (4.1). As stated in Appendix A, we assume that up to z = 1 
photometrically selected clusters are also spectroscopically confirmed, and consequently 
we assume a redshift error cTz = 0.001(1 + z) for 0.2 < z < 1, and az = 0.03(1 + z) for 
1 < z < 2. We have verified that in this case 99.95% of the signal comes from z < 1, for 
both neutrino mass and dark energy parameters wq and Wa- This is due to the damping 
effect on the observed power spectrum caused by the redshift error cr^ which increases 
linearly with redshift (see Eq. (A. 5)). Supposing instead that photometrically selected 
galaxy clusters could be spectroscopically confirmed up to z = 2, this would imply 96% of 
the signal coming from z < 1 as regards neutrino mass constraints, while for errors on the 
dark energy equation of state, one would get the 99% of the signal in terms figure of merit 
from z < 1. After 2019, future galaxy surveys will probably be able to reach spectroscopic 
precision on az even for redshifts z > 1. 
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Figure 2: 2-parameter projected 68% C.L. contours in the M^-qa subspace with = 
Us, wq, Wa, \og{10^'^ As) . The sohd black hne and the dotted blue line correspond to BAO+Pc(fc)- 
shape cluster data from a Euclid-\ike survey in combination with Planck priors for a general cos- 
mology, and a ACDM Universe, respectively. The dashed red line, and the dot-dashed orange line 
are obtained with the further addition of cluster counts data, and represent the confidence contours 
obtained for the two cosmological models, respectively. 



As we can observe, the marginalised errors of all the cosmological parameters are 1- 
2 orders of magnitude larger than the corresponding errors obtained from galaxy power 
spectrum measurements [32]. This is expected and due to the fact that the galaxy cluster 
spatial density is much lower than the density of spectroscopically selected galaxies from 
a Euclid-like survey, and evidently this effect is not compensated by the larger value of 
the effective bias. Let us notice that, having marginalised over the cluster power spectrum 
amplitude and redshift-space distortions, the information from BAO-|-Pc(fc)-shape does 
not provide any constraint on log(10^''^s). In this case the 1-a error on the total neutrino 
mass is quite large, a{M^) = 2.6 eV and, except for wq, all the cosmological parameters are 
non- negligibly correlated with M^/, with the main correlations given by — rirn, Afj, — il.},, 
and — h. 

Looking at the 4*^ and 5*^ columns of the upper panel of Table 1, we can observe that 
all the constraints improve considerably, by 1 — 2 orders of magnitude, when Planck priors 
are added to BAO-)-Pc(^)-shape cluster data, except for the dark energy equation of state, 
owing to the CMB weakness in constraining wq and Wa simultaneously. Focusing on the 
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neutrino mass, in this case (j{My) = 0.23 eV and all the parameter degeneracies with Mi, 
are broken, except for log(lO^^^s), which is now constrained by CMB measurements and 
is highly degenerate with r = 0.9, since these two parameters produce opposite effects 
on the cluster power spectrum. 

Let us now consider a less general cosmology, i.e. the ACDM model, where the cur- 
vature and the dark energy equation of state are kept fixed to their fiducial values. 
Looking at the 2"^*^ column of the lower panel in Table 1, we note that the 1-a error on the 
neutrino mass decreases by a factor of 1.5 with respect to the general case, and the com- 
bination with Planck priors (4*"^ column) gives a{My) = 0.08 eV, which is comparable with 
the lowest value of the total neutrino mass, = 0.05 eV, predicted by neutrino oscillation 
experiments. This means that the combination of CMB data with BAO-|-Pc(fc)-shape mea- 
surements from galaxy clusters, selected from future nearly all-sky galaxy catalogs, could 
detect values of the neutrino mass close to the minimum one assuming a ACDM Universe. 
Anyway, all the cosmological parameters keep a large correlation with M^, implying that 
information from external data, as e.g. a Gaussian prior on h, could help to further improve 
the neutrino mass constraint. 

5.2 COUNTS-I-CMB 

The cosmological parameter forecasts obtained from cluster counts in a Euclid-like galaxy 
survey are reported in the 6**^ and 7**^ columns of the upper panel in Table 1 for a general 
cosmology, and in the corresponding columns of the lower panel for the ACDM case. The 
four nuisance parameters described in § 2.1 have been marginalised over. The constraining 
power of cluster counts depends on the considered parameter: for example the 1-a errors 
on and wq are respectively ~ 2.5 and ~ 4 times lower than the ones obtained from 
cluster BAO-|-Pc(A;)-shape measurements, for both the model cosmologies. Nonetheless, 
the constraints on all the remaining cosmological parameters are much worse than the 
BAO+Pc(^)-shapc forecasts, and in particular the error on the total neutrino mass from 
cluster counts for the general cosmology, a{My) = 28.5 eV, is ~ 11 times larger than 
from the cluster power spectrum, while, for the ACDM model it is ~ 2.5 times larger, 
cf{Mi,) = 4.3 eV. For cluster counts as well, all the cosmological parameters are highly 
correlated with and the larger correlation coefficients r are given by Mv — Vtm, My — Q,},, 
My -h, My- Hs, and My - log(10^°As). 

Adding cluster counts data to Planck priors breaks degeneracies in CMB measure- 
ments, especially for Qk, wq and Wa, so that cosmological parameter constraints coming 
from their combination improve considerably when compared to cluster counts alone. Any- 
way, also in this case, all the cosmological parameters remain strongly correlated with the 
neutrino mass. For COUNTS-I-CMB we find a{My) = 0.26 eV for a cosmological model 
which includes curvature and an evolving dark energy component, and a{My) = 0.17 
eV for a ACDM cosmology. These findings have to be compared with the correspond- 
ing errors from cluster BAO+Pc(^)-shape measurements, given by a{My) = 0.23 eV and 
a{My) = 0.08 eV, respectively. 
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5.3 BAO+COUNTS 



Being complementary probes, the separate constraining powers of cluster counts and clus- 
ter BAO+Pc(fc)-shape data are notably amplified when combined together. In the case 
of the present work, the degeneracy breaking, provided by such a combination, is mainly 
due to the inclusion, in cluster count measurements, of information from the power spec- 
trum amplitude and from the growth of structure (encapsulated in n{M, z)), which, on the 
contrary, are marginalised over for BAO-|-Pc(^)-shape measurements. This additional in- 
formation helps to mitigate the degeneracies between the cosmological parameters, present 
when each probe is taken individually. Nonetheless, it is worth noting that dN/dz and 
Pc{k,z) describe different phenomena, abundance and clustering respectively, and have a 
different dependence on cosmological parameters; for instance, the growth rate of perturba- 
tions fg influences mainly the redshift evolution of cluster abundance, while free-streaming 
of neutrinos suppresses fluctuation power on small scales, mainly affecting the cluster corre- 
lation. Therefore, the inclusion even of amplitude and growth information in cluster power 
spectrum measurements would improve the cosmological parameter errors, but, anyway, 
would not replace the role of cluster counts in statistical constraints [76]. Moreover, these 
two probes are affected by totally different systematics, and this contributes in breaking 
the degeneracies present in each method. 

The qualitative explanation above can be verified by comparing the 10**^ column of the 
upper panel of Table 1 with the 2"*^ and 6**^ columns. We notice how the l~a errors for 
Qm, ^K, ^fe, h, and log(lO^'^^s) are reduced by ~ 1 order of magnitude with respect to the 
minimum of the corresponding errors from cluster counts and cluster BAO+Pc(fc)-shape 
measurements taken separately. The decrease is less pronounced but still important for Ug, 
Wo, Wa- In particular for the total neutrino mass we find a{M,y) = 0.9 eV in a general model 
cosmology, which is comparable to constrains from the Lya forest at 95% C.L. [24], and 
a{Mi,) = 0.57 cV in a ACDM cosmology (lower panel of Table 1), comparable to constraints 
from redshift space distortions [43] . The cosmological parameters more strongly correlated 
with are in this case Q.^, h, and log{W^^As). 

5.4 BAO+COUNTS+CMB 

In this Section we consider the total combination of cluster counts and BAO+Pc(fc)-shape 
measurements with CMB data. The results are shown in the last two columns of Ta- 
ble 1, while Figs. 1-2 represent the 2-parameter projected 68% C.L. contours in the M^-Qa 
subspace. The solid black line and the dotted blue line correspond to BAO-|-Pc(^)-shape 
cluster data from a Euclid-Yike survey in combination with Planck priors for a general 
cosmology, and a ACDM Universe, respectively. The dashed red line, and the dot-dashed 
orange line are obtained from the further addition of cluster counts data, and represent 
respectively the confidence contours obtained for the two cosmological models. These plots 
clearly display the contribution of cluster counts in constraining cosmological parameters, 
and show that the main effect results from breaking degeneracies when the cosmological 
model includes curvature and an evolving dark energy component (compare the solid black 
line with the red dashed one), so that the 1-a errors decrease by ~ 4 — 7 times with re- 
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spect to the combination CMB+BAO+i-*c(^)-shape for all the parameters, except for Ug, 
log(lO^'^^s)) and M^,. Specifically, for the total neutrino mass error in a general cosmology 
we find a{M^) = 0.18 eV, which reduces to cr(M^) = 0.076 eV in a ACDM Universe^. 
It is worth noting that, in the latter case, all the cosmological parameter constraints are 
comparable to the corresponding ones as obtained in a general cosmology from the galaxy 
CMB+BAO+Pg(A;)-shape measurements in a Euclid-like spectroscopic galaxy survey [32]. 
This means that, even if the galaxy clustering constraining power is undoubtedly superior 
to the galaxy cluster one, the two probes are complementary and their combination could 
greatly help to improve constraints on all the cosmological parameters, including neutrino 
mass and dark energy. 

Finally, the 2"*^ and 4**^ columns of Table 2 show that, fixing the four nuisance pa- 
rameters A, B, CThiMfii Siiid /3, the 1-a errors on the total neutrino mass reduce from 
0.18 eV to 0.11 eV for a general cosmology, and from 0.076 eV to 0.03 eV for a ACDM 
model. Therefore, if the true underlying cluster masses were known without uncertainties, 
it would be possible to detect at 3-a level and measure the mass of cosmic neutrinos in a 
cosmology-independent way, if the sum of neutrino masses is above 0.3 eV, while assuming 
spatial flatness and cosmological constant otherwise. 

6. Conclusions 

In this work we adopted a Fisher matrix approach to forecast constraints on a wide array of 
cosmological parameters from the photometric galaxy cluster catalog produced by future 
Euclid-like surveys. Specifically, we focused attention on the total mass of neutrinos as 
derived by a combination of cluster number counts and the shape of the cluster power 
spectrum, including BAO information. In addition to the cosmological parameters of the 
standard ACDM model and the total neutrino mass, we also considered non-vanishing 
curvature, as well as a dynamical evolution of the dark energy equation of state. We 
marginalised our constraints over the growth function, the amplitude of the scalar curvature 
power spectrum, and a set of nuisance parameters intended to represent the uncertainties in 
assigning true masses to galaxy clusters. In order to make our analysis more complete, we 
also added priors on cosmological parameters coming from the ongoing CMB experiment 
Planck. Our main results can be summarised as follows. 

• Constraints coming from cluster counts and power spectrum shape+BAO separately 
are in general relatively weak. Errors on parameter estimated using the cluster corre- 
lation function alone are up to two orders of magnitude larger than what is obtained 
with the galaxy correlation function, due to the lower spatial number density of clus- 
ters. Number counts alone do better than this only for the matter density parameter 
and the present-day dark energy equation of state, while they perform worst in all 

*For completeness and comparison reasons, we have considered also the case where only photometric 
redshifts are available. Under this assumption, for the combination CMB+BAO+Pc(fe)-shape we find 
a{Mu) = 0.24 eV in a general cosmology, which reduces to a{Mv) = 0.16 eV in a ACDM Universe, i.e. 
neutrino mass constraints would worsen by a factor ~ 1.3 and ~ 2, respectively, if a spectroscopic follow-up 
of galaxy cluster up to z = 1 would not be feasible. 
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other cases. Specifically, the error on the total neutrino mass is a{Mi,) = 2.6 eV for 
cluster correlation function alone and a{Mi,) = 28.5 eV for cluster counts alone in 
the context of a general cosmological model including curvature and evolving dark 
energy. 

• Combining the cluster power spectrum shape+BAO with cluster counts exploits the 
advantages of both approaches, and hence greatly improves the constraints of indi- 
vidual probes. The errors on cosmological parameters estimated by the combination 
of the two tests are reduced in general very substantially for all parameters, and by 
up to an order of magnitude for some of them. For a general cosmology, the error 
on the total neutrino mass is brought down to a{Mj^) = 0.9 eV, comparable with 
present constraints from the Lya forest [24] and CMB in non-flat models [16]. 

• Constraints on cosmological parameters are also greatly tightened by either i) adopt- 
ing Planck priors, and ii) reducing the number of free parameters by shrinking the 
cosmological model to a standard ACDM one. For a general cosmology, adding Planck 
priors to the cluster power spectrum shape-|-BAO reduces the parameter errors by 
1 — 2 orders of magnitude, except for the dark energy parameters, while performing 
the same operation on the cluster counts also improves the constraining power by 
more than one order of magnitude. Moreover, these priors help to remove many of 
the degeneracies between cosmological parameters, thus reducing the error on the 
total neutrino mass from either probes down to a{My) ~ 0.25 eV, which goes down 
to a{Mi,) ~ 0.08 - 0.17 if only a ACDM background is considered. 

• The combination of cluster counts, power spectrum shape+BAO, and Planck priors 
produces very competitive constraints. In a ACDM context, the errors on all cosmo- 
logical parameters are compatible with those derived for a general cosmology from 
the galaxy power spectrum shape-|-BAO performed for the spectroscopic catalog of 
a future Euclid-Yike survey in combination with Planck. Errors on the total neutrino 
mass are now down to a{Mi,) = 0.18 eV, and a{Mi,) = 0.076 cV if the cosmologi- 
cal model is enforced to be standard ACDM. The latter value is only 1.5 times the 
minimum neutrino mass admitted by neutrino oscillation experiments. 

• Fixing the value of the nuisance parameters, which equals assuming a perfect knowl- 
edge of cluster true masses, also greatly improves parameter constraints. It is quite 
unlikely that this level of precision can be reached for all clusters in the Euclid photo- 
metric catalog, while it might be a reasonable assumption for a subsample of objects 
having an extensive multi-wavelength follow-up. Nonetheless, if this will indeed be 
the case, the errors on all estimated parameters would be reduced by up to a fac- 
tor of a few. Specifically, the error on the total neutrino mass would go down to 
a{Mi,) = 0.1 eV for the general cosmology and to only a{M^) = 0.034 eV for the 
ACDM one, which means that the minimum neutrino mass could be detected. 

We emphasize that the analysis presented here is quite conservative from the point of view 
of the cosmological information, in that it does not include any guidance from the growth 



-17- 



function in tlic observed galaxy cluster power spectrum, the amplitude of the scalar cur- 
vature power spectrum, and the redshift space distortions. At the same time, we assumed 
throughout this work a Gaussian distribution of initial density and curvature perturba- 
tions. It has been repeatedly shown that the primordial non-Gaussianity has a substantial 
impact on both the number counts of massive galaxy systems and their correlation function 
[77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88]. Inclusion of this primordial non-Gaussianity 
would have the effect of adding at least another parameter to the study, hence loosening 
somewhat the constraints that have been found here [89]. 

The present investigation shows in a neat way that, whilst taken at face value the 
galaxy correlation function is superior to the cluster correlation function, the combination 
of the two, and of the latter with cluster number counts, can be of invaluable help in 
breaking the degeneracies between, and hence effectively improving the constraints on, 
cosmological parameters. Therefore galaxy clusters are expected to play a key role for 
precision cosmology with future large galaxy surveys, addressing even questions related to 
fundamental physics such as the value of the neutrino mass. 
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A. The observed galaxy cluster power spectrum 

The Fisher matrix is defined as the second derivative of the natural logarithm of the 
likelihood surface about the maximum. In the approximation that the posterior distribution 
for the parameters is a multivariate Gaussian with mean fj, = (x) and covariance matrix 
C = (xx*) — nn^, its elements are given by [59, 68, 90, 91] 



Fii = -Tr 



where x is a N-dimensional vector representing the data set, whose components Xi are 
the fluctuations in the galaxy cluster density relative to the mean in N disjoint cells that 
cover the three-dimensional survey volume in a fine grid. The {6i\ denote the cosmological 
parameters within the assumed fiducial cosmology. In the limit where the survey volume is 
much larger than the scale of any features in the observed galaxy cluster power spectrum, 
it has been shown [69] that it is possible to redefine a;„ in Eq. (A.l) to be not the density 
fluctuation in the n^^ spatial volume element, but the average power measured with the 
FKP method [92] in a thin shell of radius A;„ in Fourier space. 

In order to explore the cosmological parameter constraints from a given galaxy cluster 
survey, we need to specify the measurement uncertainties of the cluster power spectrum. 
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Therefore, the statistical error on the measurement of the galaxy cluster power spectrum 
Pc{k) at a given wave-number bin is [92] 



APr 



2(27r)^ 

Purvey ^^^^^M 



1 + ^ , (A.2) 

where Uc is the mean number density of galaxy clusters given by Eq. (2.7), ^survey is the 
comoving survey volume of the galaxy cluster survey, and /x is the cosine of the angle 
between k and the line-of-sight direction ^ = k ■ f/k. 

Actually, the observed galaxy cluster power spectrum can be different from the true 
spectrum, and it can be reconstructed assuming a reference cosmology (which we consider 
to be our fiducial cosmology) as (e.g. [60]) 



Pohs{kref±, krei\\, z) = 



Da{z)1,H{z) 
Da{zYH{zU 



Pc{kTei±, KefW , Z) + Pshot , 



(A.3) 



where 



Pc(^ref±,feref||,^) = ^e(^) 



l + /3{z,k) 



k r\\ 

ret 



'^ref± 



xPL{k,z) 



(A.4) 



In Eq. (A.3), H{z) and Da{z) are the Hubble parameter and the angular diameter distance, 
respectively, and the prefactor {DA{z)'^^fH{z))/{DA{z)'^H{z)ref) encapsulates the geometri- 
cal distortions due to the Alcock-Paczynski effect [60, 93, 94]. Their values in the reference 
cosmology are distinguished by the subscript 'rcf, while those in the true cosmology have 
no subscript. k± and /cy are the wave-numbers across and along the line of sight in the 
true cosmology, and they are related to the wave-numbers calculated assuming the refer- 
ence cosmology by krei± = k^DA{z) / DA{z\ei and fcrefn = kiiH{z\ei/H{z). Pghot is the 
unknown white shot noise that remains even after the conventional shot noise of inverse 
number density has been subtracted [60], and which could arise from clustering bias even 
on large scales due to local bias [96]. In Eq. (A.4), bc{z) is the effective bias of Eq. (2.3) 
between galaxy clusters and matter density distributions, and /3{z,k) = fg{z, k)/be{z) is 
the linear redshift-space distortion parameter [97], which in the presence of massive neutri- 
nos depends on both redshift and wave-numbers, since in this case the linear growth rate 
fg{z,k) is scale dependent even at the linear level. We estimate fg{z,k) using the fitting 
formula of Rcf. [57]. For the linear matter power spectrum Pi^{k,z), we can encapsulate 
the effect of massive neutrino free-streaming into a redshift dependent total matter linear 
transfer function T{k,z) [98, 99, 100], so that Pi,{k,z) in Eq. (A.3) takes the form 



PL{k,z) 



Sn^c'^koA, 



T\k,z) 



G{z) 



G{z = 0), 



(A.5) 



where G{z) is the usual scale independent linear growth factor in the absence of massive 
neutrino free-streaming, i.e. for — t- (see Eq. (25) in Ref. [100]), whose fiducial value in 
each redshift bin is computed through numerical integration of the differential equations 
governing the growth of linear perturbations in the presence of dark energy [101]. The 
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redshift- dependent linear transfer function T{k, z) depends on matter, baryon and massive 
neutrino densities (neglecting dark energy at early times), and is computed in each redshift 
bin using CAMB^ [46]. 

7 2 2 2 

In Eq. (A.5) we have added the damping factor e~'^ '^^ , due to redshift uncertainties, 
where = {dr/dz)az^ r{z) being the comoving distance [60, 102]. Here we adopt = 
0.001(1 + z), for 0.2<z< 1, and = 0.03(1 + z) for 1 < z < 2 [33], since we assume that 
photometrically selected clusters will be confirmed also spectroscopically for low redshifts. 

The power spectrum of primordial curvature perturbations, P-ji^k), is 



where = 0.002/Mpc, Ag = 2.525 x 10 ^ is the dimensionless amplitude of the primordial 
curvature perturbations evaluated at a pivot scale A;o, and Ug is the scalar spectral index 



Finally, to minimize nonlinear effects, in our Fisher matrix analysis we restrict wave- 
numbers to the quasi-linear regime, so that, for z < 1, femax is given by requiring that the 
variance of matter fluctuations in a sphere of radius R is cf^{R) = 0.25 for i? = tt/ (2A;max)- 
This gives /Cmax — 0.1/t Mpc"-^ at 2; = and fcmax — 0.2/i Mpc~^ at z = 1, well within the 
quasi-linear regime. In addition, for 1 < z < 2, we choose fcmax = 0.03/i Mpc~^. Finally, 
in all the calculations we impose fcmin = 10~^/i/Mpc. 

B. The Planck Fisher matrix 

In this work we use the Planck mission parameter constraints as CMB priors, by estimating 
the cosmological parameter errors via measurements of the temperature and polarization 
power spectra. As CMB anisotropies, with the exception of the integrated Sachs- Wolfe 
effect, are not able to constrain the equation of state of dark energy {wQ,Wa)^^ , we follow 
the prescription laid out by DETF [103, 104]. 

We do not include any B-mode in our forecasts and assume no tensor mode contribu- 
tion to the power spectra. We use the 100 GHz, 143 GHz, and 217 GHz channels as science 
channels. These channels have a beam of ^fwhm = 9.5', ^fwhm = 7.1', and ^fwhm = 5', 
respectively, and sensitivities of = 2.h^iK/K, ax = 2.2/.iK/K, ax = 4.8/xiC/K for 
temperature, and ap = AjiK/K^ ap = 4:.2fiK/K, ap = 9.8fiK/K for polarization, respec- 
tively. We take /sky = 0.80 as the sky fraction in order to account for galactic obstruction, 
and use a minimum ^-mode £min = 30 in order to avoid problems with polarization fore- 
grounds and not to include information from the late Integrated Sachs- Wolfe effect, which 
depends on the specific dark energy model. We discard temperature and polarization data 
at £ > 2000 to reduce sensitivity to contributions from patchy reionisation and point source 
contamination (see [103] and references therein). 

'^http://camb.info/ 

^"On the contrary, using {wo,Wa) as model parameters to compute the CMB Fisher matrix could artifi- 
cially break existing degeneracies. 




(A.6) 



[52]. 
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We assume a fiducial cosmology with an anti-corrclatcd isocurvature contribution, 
varying dark energy and curvature. Therefore we choose the following set of parame- 
ters to describe the temperature and polarization power spectra = {oJm,oJb, M,^, 100 x 
6s Aog{10^'^ As), ns,wo,Wa), where 6s is the angular size of the sound horizon at last scat- 
tering and wq and Wa are the dark energy parameters according to the CPL parametrization 
of the dark energy equation of state = wq + Wa zj (1 + z). 

The Fisher matrix for CMB power spectrum is given by [105, 106]: 



where Oi are the parameters to constrain, Cx,i is the harmonic power spectrum for the 
temperature-temperature {X = TT), temperature-E-polarization {X = TE) and the E- 

polarization-E-polarization {X = EE) power spectrum. The covariance 

COVxY of the 

errors for the various power spectra is given by the fourth moment of the distribution, 
which under Gaussian assumptions is entirely given in terms of the Cx,i with 

COVt.t = h {Ct,i + W^^B-^f (B.2) 

COY E,E = fe {Ce,i + W-'Bf'') ' (B.3) 

COYte,te = ft [C^E,1 + {Ct,i + W-'Bf^) {Ce,i + W'^Bf^) ] (B.4) 

COVt,£; = feC^E,i (B.5) 

COYt,te = feCTE,i {Ct,i + W^^Bf^) (B.6) 

COYe,te = feCTE,i {Ce,i + Wp^Bf^) , (B.7) 

where /, = 2/{{2i + l)f,^y), Wt,p = EcW^p, W^p = (a^^p^Fwhm)"' being the weight 
per solid angle for temperature and polarization respectively, with a 1-a sensitivity per 
pixel of cj^p and a beam of ^f^hm extent, for each frequency channel c. The beam window 
function is given in terms of the full width half maximum (fwhm) beam width by BJ = 
T.cm^W^.p/WT,P^ where {Bif = exp(-^(£+l)/(Z^)2), = {e^J'^ ^i8\n2) and 
/sky is the sky fraction [107]. 

We then calculate the Planck CMB Fisher matrix with the help of the publicly available 
C AMB code [46] . Finally, we transform the Planck Fisher matrix for the DETF parameter 
set to the final parameter sets q considered in this work (see § 4) , using the transformation 
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